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The complete form of the equation of state of strangeness rich proto-neutron and neutron star 
matter has been obtained. The currently obtained lower value of the Uj^ potential at the level of 5 
MeV permits the existence of additional parameter set which reproduces the weaker AA interaction. 
The effects of the strength of hyperon-hyperon interactions on the equations of state constructed for 
the chosen parameter set have been analyzed. It has been shown that replacing the strong Y — Y 
interaction model by the weak one introduces large differences in the composition of a proto-neutron 
star matter both in the strange and non-strange sectors. Also concentrations of neutrinos have been 
significantly altered in proto-neutron star interiors. The performed calculations have indicated that 
the change of the hyperon-hyperon coupling constants affects the value of the proto-neutron star 
maximum mass. 
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I. INTRODUCTION 



The mass of a star is a key factor in determining the 
stellar evolution. Low-mass stars developing a degenerate 
core can evolve to the white dwarf stage. Evolutionary 
models of massive stars, with the main sequence mass ex- 
ceeding 8M Q show that their evolution proceeds through 
sequences of consecutive stages of nuclear burning until 
the iron core is formed. The formation of 56 Fe nucleus, 
with the maximum binding energy per nucleus, indicates 
the beginning of the end of massive star's life as a nor- 
mal star. The interior composition of a massive star in 
a very advanced stage of evolution prior to core collapse 
reveal a shell structure. Each shell has different chemi- 
cal composition of gradually heavier elements. After the 
nuclear fuel is exhausted in a star, nuclear fusion in the 
central part of the star can no longer supply enough en- 



ergy to sustain a high thermal pressure. The obtained 
theoretical models of the evolved massive stars indicate 
that they develop central iron cores of mass ~ 1.5AJ© 
with the electron gas as the dominant source of pressure. 
Both thermal pressure and the pressure of degenerate 
electron gas attempt to support the core and determine 
the following equation of state of the innermost region of 
the star [h] 
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In this equation Y e = n e /ni, is the electron concentration, 
p is the density, ks is Boltzmann's constant and T is the 
adiabatic exponent. Owing to the high temperature of 
the core the adiabatic exponent is close to the critical 
value 4/3 and the core become dynamically unstable. 
The collapse of the core is accelerated by the combination 
of two mechanisms [2] : reactions of electron capture and 
photodisintegration of iron-peak nuclei into a particles 
56 Fe — ► 13a + An. Iron dissociation is an endothermic 
reaction, the energy cost equals 124.4 MeV, which comes 



from the kinetic energy of nuclei and electrons and de- 
creases their thermal pressure. The resulting pressure 
deficit is compensated by a further contraction of the 
core. Electron capture occurs on nuclei and on free pro- 
tons, diminishes the electron fraction Y e and also leads to 
pressure deficit. As a consequence the core collapse sets 
in and the star undergoes supernova explosion, ejecting 
a large part of its mass. Dynamical conditions in the 
central region of the star split the core into two separate 
parts: homologously collapsing inner core and the free 
falling outer one. After the inner core reaches the den- 
sity comparable with the nuclear density, the equation of 
state stiffens, nuclear forces are expected to resist com- 
pression and stop the collapse. The bounced inner core 
settles into a hydrostatic configuration in the dynami- 
cal time scale of milliseconds. At the boundary between 
the inner core and the supersonically falling outer core a 
pressure wave, which originates from the center, become 
a shock wave. The shock wave travels into the outer 
core through the material that is falling towards the cen- 
ter. This outgoing shock dissociates the iron nuclei into 
protons and neutrons. The shocked region attains high 
temperature at comparatively low densities. Under these 
conditions electron degeneracy is not high and relativis- 
tic positrons are thermally created. The thermal energy 
is carried out by neutrinos which are produced in the 
following reactions e + + e~ — > v + 9, e + + n — ► 9 + p. 
Within the shock and in the region behind the shock 
the velocity of the infalling matter is significantly re- 
duced and the falling matter eventually settles onto the 



2 

newly formed protoneutron star. Electron neutrinos pro- 
duced in e~ + p — > v e + n carry out electron- type lepton 
number. Dclcptonization in the outer envelope proceeds 
faster than that in the inner core. After stellar collapse 
the electron lepton number is trapped inside the matter. 
The collapse of an iron core of a massive star leads to 
the formation of a core residue, namely a proto-neutron 
star which can be considered as an intermediate stage be- 
fore the formation of a cold compact neutron star. After 
the core collapse the electron lepton number is trapped 
inside the matter. This leads to large electron neutrino 
v e chemical potential. Evolution of a proto-neutron star 
which proceeds owing to the neutrino and photon emis- 
sion causes the star change itself from a hot, bloated 
object into a cold compact neutron star and can be de- 
scribed in terms of separate phases starting from the mo- 
ment when the star becomes gravitationally decoupled 
from the expanding ejection. Through the events which 
lead up to the formation of a neutron star, two limiting 
cases can be distinguished: the very beginning stage char- 
acterized by the low entropy density s = 1 (in units of 
the Boltzmann's constant), unshocked core with trapped 
neutrinos Yl = 0.4 and low density, high entropy outer 
layer. The final stage is identified with a cold dclcp- 
tonized object (Y v = 0, s = 0). These two distinct stages 
are separated by the period of deleptonization when the 
neutrino fraction decreases from the nonzero initial value 
(Y v 7^ 0) which is established by the requirement of the 
fixed total electron lepton number at Yl = 0.4, to the 
final one characterized by Y v = 0. 



II. PROTO-NEUTRON STAR MODEL 

The presented model of a nascent neutron star is con- 
structed under the assumption that the star can be di- 
vided into two main parts: the dense core and the outer 
layer. As it was stated above, at the very beginning of 
a proto-neutron star evolution neutrinos are trapped on 
the dynamical time scale within the matter both in the 
core and in the outer layer. The estimated electron lepton 
fraction (Yl = Y e + + Y Ve ) at trapping, with the value 
of ~ 0.4 and the assumption of constant entropy, allows 
one to specify the star characteristics at the conditions 
prevailing in the star interior. Assuming the evolutionary 
scenarios presented in the paper by Prakash et al.[3] the 
following stages in the life of a proto-neutron star have 
been distinguished: 

• stage 1: the post bounce phase - the bounced inner 
core settles into a hydrostatic configuration. Under 
the assumption of very low kinetic energy of the 
matter behind the shock also the structure of the 
outer envelope can be approximated by hydrostatic 
equilibrium. Thus a proto-neutron star model im- 
mediately after the collapse can be constructed bas- 
ing on the following physical conditions: the low en- 
tropy core s = 1 with trapped neutrinos Yl = 0.4 
is surrounded by high entropy s = 2 — 5 outer layer 
also with trapped neutrinos. The energy is carried 
out from a proto-neutron star by neutrinos and an- 
tincutrinos of all flavors whereas lepton number is 
lost by the emission of electron neutrinos which are 
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produced in the j3 process e~ +p — > v e + n. Delep- 
tonization in the outer envelope proceeds faster 
than that in the inner core. 

• stage 2: deleptonization of the outer layer is com- 
pleted. The physical conditions that characterized 
the whole system are as follows: the core with the 
entropy s — 2 and Yl = 0.4 and the deleptonized 
outer envelop with s = 2 and Y v = 0. 

• stage 3: during this stage the deleptonization of the 
core takes place, after which the core has neutrino- 
free, high entropy s — 2 matter. Thermally pro- 
duced neutrino pairs of all flavors are abundant, 
and they are emitted with very similar luminosities 
from the thermal bath of the core. 

• stage 4: cold, catalyzed object. 

III. THE EQUATION OF STATE AND THE 
EQUILIBRIUM CONDITIONS 

The solution of the hydrostatic equilibrium equation 
of Tolman, Oppcnhcimer and Volkov which allows one to 
construct a theoretical model of a proto-neutron star and 
to specify and analyze its properties demands the specifi- 
cation of the equation of state. The mean field approach 
to the relativistic theory of hadrons has been used exten- 
sively in order to describe properties of nuclear matter 
and finite nuclei. The degrees of freedom relevant to this 
theory are nucleons interacting through the exchange of 
isoscalar-scalar a, isoscalar- vector w, isovector- vector p 
and the pseudo-scalar 7r mesons. In the relativistic mean 



field approach nucleons are considered as Dirac quasipar- 
ticles moving in classical meson fields. The contribution 
of the 7r meson vanishes at the mean field level. The 
chiral effective Lagrangian proposed by Furnstahl, Serot 
and Tang constructed on the basis of effective field theory 
and density functional theory for hadrons gave in the re- 
sult the extension of the standard relativistic mean field 
theory and introduced additional non-linear scalar- vector 
and vector-vector self-interactions. This Lagrangian in 
general includes all non-renormalizable couplings consis- 
tent with the underlying symmetries of QCD and can 
be considered as the one of the effective field theory of 
low energy QCD. Applying the dimensional analysis of 
Georgi and Manohar [6], [7] and the concept of natu- 
ralness one can expand the nonlinear Lagrangian and 
organize it in increasing powers of the fields and their 
I 



derivatives and truncated at given level of accuracy 
In the high density regime in neutron star interiors when 
the Fermi energy of nucleons exceeds the hyperon masses 
additional hadronic states are produced |4] , fjj] , 1 9] , [lO] , 
111]. Thus the considered model involves the full octet 
of baryons interacting through the exchange of a mesons 
which produce the medium range attraction and the ex- 
change of u> mesons responsible for the short range re- 
pulsion. The model also includes the isovector mesons p. 
In order to reproduce attractive hyperon-hyperon inter- 
action two additional hidden-strangeness mesons, which 
do not couple to nucleons, have been introduced, namely 
the scalar meson /p(975) (denoted as a* ) and the vector 
meson 0(1020) [12] . If the truncated Lagrangian includes 
terms up to the forth order it can be written in the fol- 
lowing form: 
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tub denotes baryon mass whereas m, (i = a,u), p,a* ,(f>) are masses assigned to the meson fields, M is the nucleon 
mass. There are two parameter sets presented in the original paper by Furnstahl et al |l3], |l5] Gl and G2 which have 
been determined by calculating nuclear properties such as binding energies, charge distribution and spin-orbit splitting 
for a selected set of nuclei [16). For the purposes of this paper the G2 parameter set has been chosen. Calculations 
performed with the use of this parameter set properly reproduce properties of finite nuclei and make it possible to 
compare the obtained results with the Dirac-Brueckner-Hartree-Fock (DBHF) calculations for nuclear and neutron 
matter above the saturation density. The DBHF method results in a rather soft equation of state in the vicinity of the 
saturation point and for higher densities. Calculations performed on the basis of the effective FST Lagrangian with 
the use of the G2 parameter sets predict similar, soft equation of state. Due to the fact that the expectation value 
of the p meson field is an order of magnitude smaller than that of lu meson field, the Lagrangian function ( ]3.2j ) does 
not include the quartic p meson term. In addition, as this paper deals with the problem of infinite nuclear matter 
the terms in the Lagrangian function in the paper [15] involving tensor couplings and meson field gradients have been 
excluded. The derived equations of motion constitute a set of coupled equations which have been solved in the mean 
field approximation. In this approximation meson fields are separated into classical mean field values and quantum 
fluctuations which are not included in the ground state 

a = a +So a* = a* + Sq 

As it was stated above the derivative terms are neglected and only time-like components of the vector mesons will 
survive if one assumes homogenous and isotropic infinite matter. The field equations derived from the Lagrange 
function at the mean field level arc 
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The function S(m e ff,B) is expressed with the use of the integral 
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where Jb and I3B are the spin and isospin projection of baryon B, ns is given by 
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The functions fg and are the Fermi-Dirac distribution for particles and anti-particles, respectively 

fB,B = 1 - p E[(k,rn eff , B ) T n B ]/k B T- ( 3 - 13 ) 

where fig denotes the baryon chemical potential defined as 



Pb = v B + g^BWo + g P BhBr Q + </0b/o = J k 2 + m 2 efLB + 9ujBWq + 9pBhBfo + g<pBfo (3-14) 

The obtained Dirac equation for baryons has the following form 

(ij^d^ - m e ff,B - .9ws7° w o - g<j>Bl°fo)^B = (3.15) 

with m e ff } B being the effective baryon mass generated by the baryon and scalar fields interaction and defined as 

™>eff,B =m B ~ (gaBSQ + g^BS* ). (3.16) 

In the case of a proto-neutron star matter which includes hyperon degrees of freedom the considered parameterization 
has to be supplemented by the parameter set related to the strength of the hyperon-nucleon and hyperon-hyperon 
interactions. The scalar meson coupling to hyperons can be calculated from the potential depth of a hyperon in the 
saturated nuclear matter (17 ) 

Uy = -ga-yso + g^YUJo (3-17) 

The considered model does not include S hyperons due to the remaining uncertainty of the form of their potential 
in nuclear matter at saturation density |18], [l9], [20] , |21|. In the scalar sector the scalar coupling of the A and 3 



hyperons requires constraining in order to reproduce the estimated values of the potential felt by a single A and a 
single S in saturated nuclear matter (p ~ 2.5 10 14 g/cm 3 ) 

U { A N) (po) = gaAs {po)~gu ) AW {po)-27-28MeV (3.18) 
ui N \p ) = g*~s {po) - g^~w { Po ) ~ 18 - 20MeV. 
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Table I: Chosen parameter sets. 



m CT [MeV] g a g u g p rji r] 2 r] p Co 
G2 520.254 10.4957 12.7624 9.48346 0.64992 0.10975 0.3901 2.6416 



Table II: Strange scalar sector parameters. 





9aA 


9aE 


9a*A 


9**e 


G2 weak 


6.41007 


3.33730 


3.89005 


11.64347 


strong 


6.41007 


3.33730 


7.93120 


12.45836 



Assuming the SU(6) symmetry for the vector coupling constants and determining the scalar coupling constants from 
the potential depths, the hyperon-meson couplings can be fixed. 

The strength of hyperon coupling to strange meson a* is restricted through the following relation [11 ] 



2Ui X) w 2U ( A X) . 



(3.19) 



which together with the estimated value of hyperon potential depths in hyperon matter provides effective constraints 
on scalar coupling constants to the a* meson. The currently obtained value of the U^' potential at the level of 5 
MeV [25| permits the existence of additional parameter set [26] which reproduces this weaker AA interaction. In the 



text this parameter set is marked as weak, whereas strong denotes the stronger AA interaction for (7 



(A) 



20 MeV 



[12 ] . The vector coupling constants for hyperons which arc determined from SU (6) symmetry constraints _24j] remain 
unchanged. They are of primary importance in determining the high density range of repulsive baryon potentials [24] 
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Since the dynamical time scale of neutron star evolution is much longer than the weak decay time scale of baryons 
(including strangeness violating processes) neutron star matter is considered as a system with conserved baryon 
number rife = X)s n s (B — n ,P> A, S - , S°) and electric charge being in chemical equilibrium with respect to weak 
decays. The chemical equilibrium is characterized by the correlations between chemical potentials which in the case 
of the process p + e~ <-> n + v e leads to the following relation 



(3.21) 



The similar weak process fi + v e «-> e + leads to 

Ma* Me M^e "I - M^ t f^asym ~l~ M&V ■ (3.22) 

Assuming that the muon neutrinos are not trapped inside the protoneutron star (/x,/ = 0) these relations involve 
three independent chemical potentials fi n , [i e and \x Vl , corresponding to baryon number, electric charge and lepton 
number conservation. In general, weak processes for baryons can be written in the following form 

B x + I «-» B 2 + vi (3.23) 

where B\ and B 2 are baryons, I and v\ denote lepton and neutrino of the corresponding flavor, respectively. Provided 
that the weak processes stated above take place in thermodynamic equilibrium the following relation between chemical 
potentials can be established 

Ms = 9s Mn - Qe B (Me - MO- (3.24) 

In this equation /j,b denotes chemical potential of baryon B with the baryon number qs and the electric charge q eB . 
Using the relation (|3.24| ) for A, £ and S hyperons the following results can be obtained 

MA = Ms° = M s ° = Mn, Ms- = Mh_ = Mn + (Me - M« e )> M P = Ms+ = Mn - (Me - M^J- ( 3 - 25 ) 

The presented constrains of charge neutrality and /3-equilibrium imply the presence of leptons which are introduced 
as free particles. Thus the Lagrangian of free leptons has to be added to the Lagrangian function Q3.2]) 

^ = E - m *)^ + E ^w^- ( 3 - 26 ) 

i 



Muons start to appear in neutron star matter in the pro- 
cess e~ «-> /i _ +^ + ^ e after // u has reached the value 
equal to the muon mass. The appearance of muons not 
only reduces the number of electrons but also affects the 
proton fraction. Numerical models of the collapsing in- 
ner core have shown that there are no muon neutrinos 
= 0) and only electron neutrinos are trapped in the 
matter [3] . After the deleptonization of a proto-neutron 



star matter there are no trapped neutrinos either and 
fi Ve = 0. In this case the number of independent chemi- 
cal potentials reduces to two (/i n and /i e ). 

IV. RESULTS 

More realistic description of a neutron star requires 
taking into consideration not only the interior region of 
a neutron star but also the remaining layers, namely the 



The pressure in the outermost layer of a proto-neutron 
star is determined by nucleons, electrons and neutrinos. 
The theoretical model of the dense inner core can be con- 
structed under the following assumptions: the matter in- 
cludes the full octet of baryons interacting through the 
exchange of meson fields, the composition is determined 
by the requirements of charge neutrality and generalized 
(3 equilibrium. The loss of lepton number from the col- 
lapsed star proceeds in separated stages. First, very fast, 
carried on in a very short time in comparison with the 
Kelvin-Helmholtz neutrino cooling time, the deleptoniza- 
tion of the surface layer takes place [3j. After that the 
deleptonization of the core proceeds by the emission of 
electron neutrinos which are produced efficiently via the 
/3-process 



inner and outer crust and surface layers. The composite 
equation of state has been constructed by adding Bonn 
and Negele-Vautherin (NV) equations of state (describ- 
ing the inner crust) to the one determined with the use 
of the G2 parameterization of the FST model ([l3],|14]) 
allows to calculate the neutron star structure for the en- 
tire neutron star density span. The proto-neutron star 
models have been constructed on the assumption that 
they include a core and an outer layer. The outer re- 
gion of a proto-neutron star is characterized by the fol- 
lowing conditions: it is less dense and less massive than 
the core, the assumptions of nearly complete disintegra- 
tion of nuclei into free nucleons and (3 equilibrium have 
been made. Chemical potentials of the constituents of 
the matter are evaluated by the requirement of charge 
neutrality and (3 equilibrium. At sufficiently low density, 
in the absence of interactions which alter the nucleon ef- 
fective mass nucleons are considered as non-relativistic. 
I 

The events between the two distinguished moments depend on details of the considered models, especially on the 
equation of state. The initial moment t = characterizes the stage with neutrinos trapped (Yj, = 0.4) both in the 
core and in the outer envelope. Introducing a parameter a, which determines the degree of deleptonization of the 
proto-neutron star matter, the moment t = can be identified with a = 1 whereas the final completely deleptonized 
stage (Yy = 0) is described by a = 0. The density at the core-outer layer interface is time dependent. This dependance 
can be presented by comparison of the interface between the inner core equation of state and the equations of state 
which have been obtained for different values of the parameter a in the outer layer. Results are shown in Fig. CD The 
deleptonization leads to the substantial reduction in the extension of the proto-neutron star envelope. The obtained 
form of the equation of state for the strong and weak Y — Y interactions are presented in Fig. [21 Both models 
have yielded identical results up to a certain value of density: the hot, neutrino-trapped matter is described by the 
stiffest equation of state whereas the hot, deleptonized and cold (T = 0) proto-neutron star models result in softer 
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Figure 1: The density dependance of the interface between the inner core equation of state and the equations of state which 
have been obtained for different values of the parameter a in the outer layer for the the strong (left panel) and weak (right 
panel) Y — Y interaction models. 




Figure 2: The equation of state of the hyperon star matter for the G2 parameter set. In the left and right panels the results 
for the strong and weak Y — Y interactions are presented. 

equations of state. However, the strong Y — Y interaction model gives for the hot, deleptonized and for the cold T = 
proto-neutron star matter equations of state stiffer than that for the weak model. 

The composition of hyperon star matter as well as the the threshold density for hyperons are altered when the strength 
of the hyperon-hyperon interaction is changed. Fig. [3] presents fractions of particular strange baryon species Yg as 
a function of baryon number density rib for the weak and strong model of the G2 parameterization, for the three 
selected evolutionary phases. All calculations have been done on the assumption that the repulsive X interaction 
shifts the onset points of £ hyperons to very high densities and they do not appear in the considered proto-neutron 
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Figure 3: Relative concentrations of A and H hyperons in hyperon star matter for the G2 parameter set. In the left and right 
panels the results for the strong and weak Y-Y interaction are presented. 



and neutron star models. The onset of S° hyperons takes place at high densities and they are not presented in 
Fig. [31 In the case of strong and weak Y — Y interactions the onset of A hyperons is followed by the onset of S~ 
hyperons. In general, when A hyperons appear they replace protons and, in the consequence, lower the energy of the 
system. This fact, through the requirement of charge neutrality of neutron star matter, results in diminished electron 
chemical potential. The populations of A hyperons obtained in the weak models in the considered evolutionary stages 
are reduced in comparison with those calculated for the strong Y — Y interaction models. Thus the weak model 
with the reduced A hyperon population leads to the increased value of electron and proton concentrations. Results 
are presented in Fig. \4\ where relative concentrations of electrons and protons are depicted. One can also compare 
the concentrations of neutrinos in the hot, neutrino trapped matter. In Fig. [5] neutrino concentrations for the two 
different models are shown. The weak model leads to lower concentrations of neutrinos. For comparison results for 
non-strange proto-neutron star matter is included. 

Once the equations of state have been calculated for each evolutionary phase and for the strong and weak models 
the corresponding hydrostatic models of proto-neutron stars have been obtained. Solutions of the Oppenheimer- 
Tolman-Volkov equation for the considered parameter sets are presented in Fig. [6} In agreement with the introduced 
evolutionary stages of a proto-neutron star the obtained mass-radius relations include a sequence of models calculated 
for a given equation of state characteristic for a specific phase of a proto-neutron star evolution. In the left and right 
panel of Fig. \6\ results for six equations of state have been plotted. The four mass-radius relations are constructed 
for the strangeness-rich baryonic matter and they represent the evolution of a proto-neutron star starting from the 
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Figure 4: Relative concentrations of electrons (the left panel) and protons (the right panel) in hyperon star matter for the G2 
parameter set. 
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Figure 5: Left panel: Density dependance of neutrino concentrations. Right panel: Neutrino concentrations as a function of 
stellar radius R. 



moment when the proto-neutron star can be modelled as a low entropy core surrounded by a high entropy outer layer. 
Neutrinos are trapped both in the core and in the outer layer. This very beginning phase of a proto-neutron star is 
followed by two subsequent stages connected with the process of dclcptonization of a proto-neutron star. The first 
during which the deleptonization of the outer layer takes place whereas in the core neutrinos are still trapped and the 
second one representing a hot deleptonized object with thermally produced neutrino pairs of all flavor abundant in 
the core and in the outer layer. The final case is exemplified by a solution obtained for cold neutron star matter which 
includes hyperons. For comparison the mass-radius relations for two limiting cases of proto-neutron star evolution 
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Figure 6: The mass-radius relations for subsequent stages of the proto-neutron star evolution. The left panel presents results for 
G2 parameterization with the strong Y — Y interaction, the right panel for the weak interaction model. Dotted lines represent 
solutions obtained for non-strange baryonic matter. 



Table III: Proto-neutron and neutron star parameters obtained for the maximum mass configurations for the G2 parameteri- 
zation with the strong Y — Y interactions. 





p/po 


R (km) 


M(M ) 


M B (M ) 


T=0 


A 


5.00 


12.73 


1.298 


1.407 


B 


21.76 


8.22 


1.320 


1.634 


s=2 


8.76 


13.528 


1.670 


1.859 



for nonstrange baryonic matter have been presented. A similar sequence of the mass-radius relations have been 
constructed for the G2 parameterization supplemented by the weak Y — Y parameter set for the strange sector. There 
are qualitative changes in proto-neutron star parameters which occur for a given evolutionary stage. The G2-strong 
parameterization leads to proto-neutron neutron stars with the reduced value of the maximum mass. For the strong 
hyperon-hyperon interaction strength in the case of the cold neutron star model besides the ordinary neutron star 
branch there exists additional stable branch of solutions which are characterized by a similar value of the mass but 
with a significantly reduced radius. For the purpose of this paper A denotes the maximum mass configuration of 
the ordinary neutron star branch and B the additional maximum. Proto-neutron and neutron star parameters are 
summarized in Tables IIII and |IVl The presented results have been obtained for maximum mass configurations. 
The analysis of the existence and stability of the additional branch in the mass-radius relation will be the subject of 
further investigations. 

Neutron stars are purely gravitationally bound object. The gravitational binding energy of a relativistic star is defined 
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Table IV: Proto-neutron and neutron star parameters obtained for the maximum mass configurations for the G2 parameteri- 
zation with the weak Y — Y interactions. 





p/po 


R (km) 


M (M ) 


M B (M ) 


T=0 


8.00 


11.63 


1.527 


1.749 


s=2 


9.12 


13.14 


1.754 


1.916 




Figure 7: The neutron star mass versus the baryon mass Mb = Mc 2 Nb- In the left and right panels the results for the strong 
and weak Y-Y interaction, respectively, are presented. 

as the difference between its gravitational and baryonic masses. 

E b<g = (M - M B )c 2 (4.2) 

where 

M = 4vr C ' drr 2 (l - 2G ^ )~h p ( r ) (4.3) 
Jo cr 

and Mb — n\,m,B is given by the total baryon number n^. The numerical solution of the above equation has been found 
for the selected equations of state for the strong and weak interactions. The results are shown in Fig. [7] which also 
includes for comparison solutions for non-strange proto-neutron and neutron star models. The obtained mass-radius 
relations for the non-strange baryonic matter confirm the well-known fact that in this case the maximum mass of a 
proto-neutron star with trapped neutrinos is lower than that of cold deleptonized matter. This has a consequence 
for the possibility of a black hole formation during the relatively long lasting phase of deleptonization. In the case of 
non-strange matter a proto-neutron star is not able to achieve the unstable configuration due to deleptonization. In 
general, neutrino trapping increases the value of the maximum mass for the strangeness rich baryonic matter. In all 
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evolutionary stages there are configurations which due to deleptonization go to the unstable branch of proto-neutron 
and neutron stars. 

The gained solutions of the structure equations allows as to carry out a similar analysis of the onset point, abundance 
and distributions of the individual baryon and of lepton species but now as functions of the star radius R. The 
analysis includes results obtained for two extreme evolutionary phases: the very beginning hot, neutrino-trapped 
matter stage and the cold, deleptonized one. Two characteristic configurations have been considered for the T = 
solution. Namely, the one connected with the maximum mass configurations marked as A and the other for the B 
maximum mass configuration of the G2 strong parameterization. The compact hyperon core which emerges in the 
interior of the maximum mass configuration consist of A, S~ and S° hyperons. The hyperon population is reduced 
to A and S~ in the case of hot neutrino trapped matter for the weak model. The process of deleptonization and 
cooling diminishes the concentration of A hyperons and increases concentrations of H _ , S° hyperons. The reduction 
of A hyperon population in the weak model is larger then that obtained for the strong Y — Y interaction model. The 
innermost hyperon core has larger radius and contains much more negatively charged 5 hyperons in the weak model. 
The relative fractions of hyperons in the core of the maximum mass configurations are presented in Fig. \S\ and Fig. 
[9} Horizontal lines in the left panels of Fig. [8] and Fig. \9\ correspond to the threshold density of individual hyperons 
in the considered model. The configuration marked as A does not contain hyperons. The appearance of S - hyperons 
through the condition of charge neutrality affects the lepton fraction and causes a drop in their contents. Charge 
neutrality tends to be guaranteed with the reduced lepton contribution. Results are shown in FigJTOl The relative 
hadron-lepton composition in this model can be also analyzed through the density dependence of the asymmetry 
parameter fa which describes the neutron excess in the system and the parameter fs — (Na + 2N-£- + 2N^o)/nb 
which is connected with the strangeness contents. FigJTl] and [12] present both parameters as functions of baryon 
number density and the stellar radius R. As the density increases, the asymmetry of the matter decreases and the 
parameter fs increases for all the considered cases. Deleptonization and cooling leads to stellar matter which is more 
asymmetric and possesses substantially enhanced strangene ss contents. 



V. CONCLUSIONS 

In this paper the complete form of the equation of 
state of strangeness rich proto-neutron and neutron star 
matter has been obtained. The considered models are 
constructed on the assumption that the proto-neutron 



star consists of two main parts: the core and the outer 
layer. The pressure in the outer envelope is determined 
by free nucleons, electrons and neutrinos in [3 equilib- 
rium. The theoretical model of the dense inner core 
is described by the Lagrangian which includes the full 
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Figure 8: The fraction of species i, Yi in the maximum mass configuration as a function of star radius for the G2 strong 
parameter set (fci? means the effective Fermi momentum which specifies the value of in equation (|3. 14| ) . 



octet of baryons interacting through the exchange of me- 
son fields. The chosen parameter set based on the ef- 
fective field theory includes nonlinear scalar-vector and 
vector- vector interaction terms which leads to rather soft 
equation of state. This is in agrement with the DBHF 
results. The meson sector has been extended not only 
by nonlinear terms but also by two additional hidden- 
strangeness mesons which reproduce attractive hyperon- 
hyperon interactions. Thus the chosen parameter set 
G2 has been supplemented by the parameter set related 
to the strength of the hyperon-nucleon and hyperon- 
hyperon interactions. The currently obtained lower value 
of the uj^ potential at the level of 5 MeV permits the ex- 
istence of additional parameter set which reproduces this 



weaker AA interaction. The main goal of this paper was 
to study the influence of the strength of hypcron-hyperon 
interactions on the properties of the proto-neutron star 
matter and through this on a proto-neutron star struc- 
ture during selected phases of its evolution. The presence 
of hyperons in general leads to the softening of the equa- 
tion of state. The behavior of the equation of state is 
directly connected with the value of the maximum star 
mass. Equilibrium conditions, namely charge neutrality 
and /^-equilibrium, determine the composition of the star. 
It has been shown that replacing the strong Y — Y inter- 
action model by the weak one introduces large differences 
in the composition of a proto-neutron star matter both 
in the strange and non-strange sectors. There is a con- 
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Figure 9: The fraction of species i, Yt in the maximum mass configuration as a function of star radius for the G2 weak parameter 
set. 




02468 10 02468 10 



R ( km ) R ( km ) 



Figure 10: Relative concentrations of A and H hyperons in hyperon star matter for the G2 parameter set. In the left and right 
panels the results for the strong and weak Y-Y interaction, respectively, are presented. 



siderable reduction of A hyperon concentration whereas 
the concentrations of S~ and S° hyperons are enhanced 
during the deleptonization. In addition, the population 



of A hyperons obtained in the weak models in the con- 
sidered evolutionary stages is lower in comparison with 
those calculated for the strong Y — Y interaction models. 
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n b ( fm " ) R ( km ) 



Figure 11: The parameter fa as a function of baryon number density for the strange and non-strange baryonic matter calculated 
for each specified evolutionary phases. The hot models with trapped neutrinos are marked as s = 2, whereas hot, deleptonized 
models as Y v — 0. In the right panel the parameter fa as a function of stellar radius for the weak and strong models are 
presented. 




n b ( fm " ) R(km) 



Figure 12: The parameter fa as a function of baryon number density for the strange and non-strange baryonic matter calculated 
for each specified evolutionary phases. The hot models with trapped neutrinos are marked as s = 2, whereas hot, deleptonized 
models as Y v — 0. In the right panel the parameter fa as a function of stellar radius for the weak and strong models are 
presented. 



Thus the weak model with the reduced A hyperon pop- 
ulation permits larger fractions of protons and electrons 
and leads to lower concentrations of neutrinos. The G2- 
strong parameterization leads to proto-neutron neutron 
stars with the reduced value of the maximum mass. For 



the strong hyperon-hyperon interaction strength in the 
case of the cold neutron star model besides the ordinary 
neutron star branch there exists additional stable branch 
of solutions which are characterized by a similar value 
of the mass but with significantly reduced radius. This 



may be connected with the _existence of the third fam- 
ily of stable compact stars [28j. The analysis of baryon 
and lepton concentrations as a function of stellar radius 
R shows that for the chosen parameter set there is no 
strange baryons in the maximum mass configuration A. 
This star resembles an ordinary neutron star. The ad- 
ditional stable configurations of neutron stars which ap- 
pear on the mass-radius diagram include large fraction of 
hyperons (hyperstar). Analysis of the concentrations of 
leptons (electrons and muons) in the compact star B as 
a stellar radius R leads to the conclusion in the leptons 
are gathered? in the outer layer of the star. Thus, the 
star has a lepton rich outer part and almost completely 
deleptonized inner core. This analysis shows that there 
exists a very strong correlation between the value of the 
maximum neutron star mass and the strength of hyperon 
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coupling constants. The inclusion of additional nonlinear 
meson interaction terms which modify the high density 
behavior of the equation of state together with the strong 
hyperon-hyperon interaction lead to the existence of ad- 
ditional stable stellar configurations with similar masses 
and smaller radii than an ordinary neutron star. The re- 
duction in radius is of the order of 2.5 km. The internal 
composition of this additional neutron star configurations 
is almost completely free of leptons. Transmutation simi- 
lar to the phase transition from the ordinary neutron star 
to the more compact hyperstar may be the main origin 
of the short gamma ray burnst [29] . 
Employing the data concerning the estimated value of 
the A well depth Uaa — 5 MeV the existence of the addi- 
tional stable branch of very compact stars have not been 
confirmed. 
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